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Abstract. In recent years there has been considerable progress in our understanding of the nature and properties of 
the reionization process. In particular, the numerical simulations of this epoch have made a qualitative leap forward, 
reaching sufficiently large scales to derive the characteristic scales of the reionization process and thus allowing for realistic 
observational predictions. Our group has recently performed the first such large-scale radiative transfer simulations of 
reionization, run on top of state-of-the-art simulations of early structure formation. This allowed us to make the first realistic 
observational predictions about the Epoch of Reionization based on detailed radiative transfer and structure formation 
simulations. We discuss the basic features of reionization derived from our simulations and some recent results on the 
observational implications for the high-redshift Ly-0£ sources. 
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INTRODUCTION 

The reionization of the universe was the last global tran- 
sition of the Intergalactic Medium (IGM), caused by the 
radiation of the First Stars. It completely transformed the 
IGM, from fully-neutral and cold to almost fully-ionized 
and hot (Tjgm ~ 10 4 K). Reionization started with the 
formation of the first sources of ionizing radiation at very 
high redshifts z ~ 30 — 40 and finished around z ~ 6 — 7. 
Currently there are still very few direct observational 
constraints on this epoch. The lack of Gunn-Peterson 
trough in the spectra of high-redshift sources implies low 
mean neutral fraction xhi <; 10 -4 out to redshift z ~ 5.5, 
which appears to rise at higher redshifts, possibly indi- 
cating that at these redshifts we are seeing the tail-end 
of reionization. The second direct observable, the inte- 
grated electron- scattering optical depth, T es , was found 
to be 0.09 ± 0.03 based on the WMAP 3-year data 0]. 
This range of values points to an early start and extended 
period of reionization. However, the error bars remain 
fairly large and thus the optical depth by itself does not 
put very stringent constraints on reionization. 

This relative lack of observational data is set to change 
dramatically in the next few years, however. A num- 
ber of large observational projects are currently under 
development, e.g. observations at the redshifted 21 -cm 
line of hydrogen [e.g. 0, H, 0, S Sl» detection of small- 
scale CMB anisotropies due to the kinetic Sunyaev- 



Zel'dovich (kSZ) effect [e.g.MJ 
redshift Ly-a emitters [e.g. [K 



9], and surveys of high- 
Ill]. The planning and 



success of these experiments relies critically upon un- 



derstanding the large-scale topology of reionization, i.e. 
the size- and spatial distribution of the ionized and neu- 
tral patches, which is best derived by large-scale simula- 
tions. Recently we presented the first large-scale, high- 
resolution radiative transfer simulations of cosmic reion- 
ization |[12i |5|, [l3|] . Here we summarize our basic con- 
clusions about the character of the reionization process, 
along with some of our recent results on its observability. 



BASIC FEATURES OF REIONIZATION 

The key results for the nature and features of reionization 
could be summarized as follows: 

• Reionization proceeds in an inside-out fashion, 
whereby the high density regions are ionized first 
and the deepest voids - last, due to the galaxies pref- 
erentially forming in and around the high density 
peaks of the density distribution lfl2ll . 

• The process is thus highly-inhomogeneous and 
patchy on large scales, with very complex geome- 
try. It is strongly modulated by the long-wavelength 
density fluctuations. As a consequence, small- 
volume simulations, which do not include these 
large-scale fluctuations underestimate the source 
biasing and the reionization patchiness fl2[ l5[[l4ll. 

• The reionization parameters, like source efficiencies 
(dependent on the stellar IMF, star formation effi- 
ciencies and escape fractions) and small-scale gas 
clumping, are only weakly constrained at present. 
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FIGURE 1. Position-redshift slices from one of our simulations, illustrating the large-scale geometry of reionization and the 
significant local variations in reionization history as seen at the redshifted 21 -cm line. Shown is the decimal log of the differential 
brightness temperature. 




FIGURE 2. The reionization history of a high density peak. The images are centered on the most massive (at z = 6) halo in our 
computational volume and are of size 100/* Mpc to the side. The snapshots are (left to right): z = 12.9, z = 9.0, z = 8.0, and 
z = 6.6. The underlying cosmological density field (dark,green) is superimposed with the ionized fraction (light, orange) and the 
ionizing sources (dark, blue dots). 



Both the overlap epoch, z v and the integrated opti- 
cal depth T es are readily reproduced by_a wide range 
of reasonable reionization scenarios 

• The H II regions have typical scales of 5-20 co- 
moving Mpc (dependent on still-uncertain param- 
eters). Although the high-redshift sources are typi- 
cally very low-mass, their strong clustering results 
in the formation of such large ionized patches (see 
Figures [T]and [2]). These characteristic scales are di- 
rectly reflected in many of the predicted reioniza- 
tion observables SB ESQ. 

• Reionization is strongly self-regulated through 
Jeans-mass filtering of the low-mass sources: more 
efficient low-mass sources (e.g. due to top-heavy 
IMF) result in more suppression of the same 
sources, thereby largely cancelling the effect of the 
higher efficiency, and vice versa H 1 311 . 

• This process of self-regulation has the effect that 
T es and z ov , the overlap redshift are only loosely re- 
lated and do not impose very strong constraints on 
the possible reionization scenarios. The reason for 
this is that z ov depends on the abundances and effi- 
ciencies of the high-mass sources, whose formation 
is not suppressed by reionization, unlike T es , which 
depends on both types of sources. Varying the ion- 
izing efficiencies of low-mass sources yields a wide 
range of T es values for the same z ov 111311 . 



HIGH-REDSHIFT LY-a SOURCES 

Observations of the high-redshift Ly-a sources, galaxies 
and QSO's, have already provided us with a wealth of 
information about the state of the IGM and the nature of 
the luminous galaxies at the end of reionization and still 
hold a lot of promise for the future. In addition to probing 
the IGM and the source luminosity function (and thus, 
indirectly, the halo mass function and ionization source 
properties), they may also be used to constrain the reion- 
ization geometry. The first rare objects form in the high- 
est peaks of the density field. The statistics of Gaussian 
fields predicts that density peaks are strongly clustered, 
especially at high redshift. As a consequence, each high- 
redshift, massive galaxy was surrounded by numerous 
smaller ionizing sources. Self-consistent simulations of 
such regions require following a large enough volume 
while at the same time resolving all the low-mass halos 
which drive reionization, which has not been possible un- 
til recently. 

In Figure [2 we illustrate several stages of the reion- 
ization history of a high density peak. The most massive 
source in our volume (at z = 6) is shifted to the centre 
using the periodicity of the computational box. At red- 
shift z= 12.9 the H II region is small and the source 
is invisible due to the damping wing due to the neutral 
gas outside. By redshift z = 9 many more haloes have 
formed, most of them in large clustered groups. The H II 




FIGURE 3. Ly-a sources at redshift z = 6.0: (left and middle) Sample LOS vs. A/comoving distance from the most massive 
galaxy. Shown are (left) the optical depth (solid), neutral fraction x 10 4 (dotted) and density (dashed), and (right) the corresponding 
transmission. The vertical lines show the position of the central source (in redshift space, i.e. accounting for its peculiar velocity), 
while the horizontal lines indicate the optical depth equivalent to 1% transmission. Redshift- space distortions due to the local 
peculiar velocities and the Ly-a line damping wing are included; and (right) Luminosity function without (black) and with 
absorption included (red). For reference, the green, dotted line shows the result of suppressing each source by 50%, which would 
be the case if e.g. all of the blue wing of the emission line, but none of the red wing is absorbed. 



region surrounding the central peak is among the largest, 
but the central source emission is still strongly affected 
by damping. Only by redshift z = 8 is the ionized region 
large enough to render the source unaffected by damping 
and potentially visible. The reionization geometry be- 
comes quite complex, most ionized bubbles are intercon- 
nected, but large neutral patches remain between them. 
Finally, by the nominal overlap z = 6.6 all ionized re- 
gions have merged into one topologically-connected re- 
gion, although substantial neutral patches remain inter- 
spersed throughout our volume, which remains largely 
optically-thick to both Ly-a and ionizing continuum ra- 
diation. Only by z ~ 6 this volume becomes on average 
optically-thin to ionizing radiation. 

In Figure [3] (left and center) we show some sample 
spectra for the same luminous source. The spectra ex- 
hibit extended high-transmission (10-60% transmission) 
regions in the highly-ionized proximity zone of the lumi- 
nous source, within 5 Mpch~ l (~ 20 A). The center of 
the peak itself is optically-thick due to its high density. 
The infall around the central peak blue- shifts photons, 
resulting in some absorption behind the redshift- space 
position of the source. Away from the proximity region 
the absorption is largely saturated, but there are a num- 
ber of transmission gaps with up to a few per cent trans- 
mission. In Figure [3] (right) we show the Ly-a source 
luminosity function at the same redshift. For the weaker 
sources roughly half of the intrinsic luminosity is trans- 
mitted (the red wing of the line), while the most luminous 
sources suffer from additional absorption due to the gas 
infall that surrounds them. Future work will quantify the 
statistics of these features and its evolution. 
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